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ABSTRACT 

The surface brightness profile of Ha emission in galaxies is generally thought to be 
confined by a sharp truncation, sometimes speculated to coincide with a star formation 
threshold. Over the past years, observational evidence for both old and young stellar 
populations, as well as individual H II regions, has demonstrated that the outer disk 
is an actively evolving part of a galaxy. To provide constraints on the origin of the 
aforementioned Ha truncation and the relation of Ha emission in the outer disk to the 
underlying stellar population, we measure the shape of the outer Ha surface brightness 
profile of 15 isolated, edge-on late-type disk galaxies using deep, long-slit spectroscopy. 
Tracing Ha emission up to 50% beyond the optical radius, i?25, we find a composite 
Ha surface brightness profile, well described by a broken-exponential law, that drops 
more steeply in the outer disk, but which is not truncated. The stellar continuum 
and Ha surface brightness both exhibit a break at ~ 0.7i?25, but the Ha emission 
drops more steeply than the stellar continuum beyond that break. Although profiles 
with truncations or single exponential laws correctly describe the Ha surface bright- 
ness profiles of some individual galaxies, flexible broken-exponentials are required in 
most cases and are therefore the more appropriate generic description. The common 
existence of a significant second surface brightness component beyond the Ha break 
radius disfavors the hypothesis that this break is a purely stochastic effect. 

Key words: galaxies: kinematics and dynamics - galaxies: formation - galaxies: 
evolution - galaxies: spiral - galaxies: structure 



1 INTRODUCTION 

The outer disks of galaxies are increasingly being recognized 
as a diverse environment whose properties may aid us in 
understanding the for mation of galaxies. Since the radio ob- 
servations of the 1 950s (|van de Hulst. Raimond. fc WoerdenI 
I1957I : lDieteijll96^ ). we have known that the baryonic com- 
ponent of spiral galaxies, as traced by neutral hydrogen, is 
significantly larger than suggested by classical size indicators 
such as (the radius where the surface brightness drops 
to 25 mag arcsec"^; we use measurements of the B-band 
surface brightness obtained from the literature throughout 
this analysis). Optical studies, however, have focused on the 
inner parts of galaxies, where, among other aspects, the ra- 
dial distribution of Ha emission has b een subject of intense 
scruti ny for many years. For example, iMartin fc KennicuttI 
(|200ll ) used narrow-band imaging of nearby galaxies to con- 



struct Ha surface brightness profiles. Such profiles usually 
exhibit a "truncation" in the Ha surface brightness at rel- 
atively large radii, although usually still within the op- 
tical radius R25. This truncation is also reflected in the 
scarcity of Ha rotation curves that extend to radii much 
beyond ^ .7 R25 in comprehensive long-slit surveys (e.g., 
IVogt et"all (|2004 ). being one of the deepest such surveys, 
observe no rotation curves beyond R25 in a sample of 329). 
A widely accepted hypothesis for such a truncation is the 
existence of a threshold in the surface g as density, below 
which star formation beco mes inefficient (|Kennicuttl Il989l : 
iMartin fc KennicuttI l200lh . Alternate explanations include 
the hypothesis that the truncation is indicative of an actual 
break in the mass distribution that is related to the initial 
formation conditions of the disk (van dor Kruit 1987). 

However, "truncation" may be too strong a descrip- 
tive statement because it is now known that low-level 
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star formation occurs at large galactocentric radii, well be- 
yond the radius corre s pondin R to the critical t hresho ld. 
iBland-Hawthorn et all (|l997t ): iFerguson et"al] (|l998al lbh 
have detected individual Ha emission regions in the outer 
disks of spiral galaxies as far as 2 R25 , a nd UV emission ap- 
pears relatively common in o uter disks jThilker et al ] |2005l : 
IZaritskv fc ChristleinI l2007l : iThilker et all l2007h . Further- 
more, in a study of NGC 38 1 4 usin g deep two-band optical 
imagmg, iHerbert-Fort et all (j2009( ) have found statistically 
significant overdensities of marginally resolved sources in the 
outer disk that are likely to be star clusters. 

This outer-disk star formation has raised new interest 
in the hypothesis of a star formation threshold. The absence 
of evidence for a break in the surface brightness profiles of 
UV emi s sion, w hich is also a star formation indicator, led 
iBoissied (|2007al lbh to suggest that the star formation rate 
does not exhibit a break, and that the Ha truncation is a 
stochastic effect. As the expectation number of star forma- 
tion regions with stars massive enough to generate a Stroem- 
gren sphere drops below unity, the Ha emission profile goes 
to zero, while low level(mass) SF continues as measured with 
the UV. Other, perhaps more exciting, possibl e explanations 
inclu de a change in the initial mass function (|Meurer et al.l 
I2OO9I ). These arguments depend critically on reliable mea- 
sures of Ha and continuum profiles out to large radii. 

At the same time, studies of faint optical contin- 
uum emission have shown that stellar disks can also be 
traced far into the outer disk. While some galaxies ex- 
hibit a single exponenti al surface brightness profile to the 
largest measurable radii (|Bland-Hawthorn et al.|[2005l ). oth- 
ers are more accurately described by broken exponen- 
tial profiles with a characteristic break radius, beyond 
which the stellar surface brightness p rofile may be shal- 
lower or steeper th an in the inner disk jPohlen et al.ll2002l : 
lErwin. Pohlen fcB cckman 2008,). These profiles are some- 
times descri bed as exponential, sub-exponential, and s uper- 
exponential (jVlaifc. Bland-Hawthorn fc Freemanll2009l ): the 
exponential and sub-exponential types are also variou sly re- 
ferred to as Freeman Types I and II |Freemanlll970l ). and 
the designation "Type HI" has come into usage to describe 
the super-exponential (up-bending) shape. 

How these outer stellar structures are related to pos- 
sible structures in the gaseous disk, and in particular. Ha 
emission, is a crucial question in understanding the origin 
of the outer disk stellar populations. Is the truncation in 
Ha emissio n a real indicato r of a truncation in star forma- 
tion, or, as iBoissieil (|2007al lbl) suggest, merely a stochastic 
effect? If, on the other hand, the break in the star formation 
surface density is real, is it also responsible for the charac- 
teristic break in the stellar continuum profile? In that case, 
is there significant star formation beyond the break, and is 
it enough to have created the outer disk stellar content in 
situ? Or must other processes be invoked to populate the 
outer disks with stars? 

To answer these questions, a systematic and quantita- 
tive study of the distribution of Ha emission at large radii 
is required. The traditional technique of narrow-band imag- 
ing with subsequent subtraction of broad-band continuum 
emission is generally insufficiently sensitive to probe to large 
radii. The limitation lies not only in the achievable signal- 
to-noise ratio (which is limited because typical narrow-band 
filters are much wider than the Ha emission line), but also 



in the stellar continuum subtraction. At large radii, because 
the spectral flux density at the peak of the Ha line is of the 
same order of magnitude as the stellar continuum, the wave- 
length dependence of the stellar continuum spectral energy 
is sufficient to significantly degrade the subtraction. For ex- 
ample, absorption troughs, which the Ha emission line is 
often embedded in, can render it undetectable. 

The way around these problems is to use higher spec- 
tral resolution. This can involve very narrow filter band- 
passes, which is now possible over a wide range of red- 
s hifts due to the increasing availability of tunab le filters 
(jBland-Hawthorn fc Jone^ Il998l : ICepa et al.lll998l ). or tra- 
ditional spectroscop y. In a pioneering effort in th e study of 
gaseous outer disks. iBland-Hawthorn et al.l (|l997l ) detected 
Ha emission at ~ 1.25 R25, beyond the truncation radius 
of the neutral hydrogen disk, using the Fabry-Perot staring 
technique. However, true spectroscopic observations provide 
improvements in the subtraction of the stellar continuum, of 
the [Nil] emission lines, and of Ha absorption troughs, al- 
lowing one to reach sensitivity li mits fainter than 10~^* erg 
s"'^ cm~'^ arcsec"^ (see also Mad sen et al.l fSooJ)). 

We have observed late-type, edge-on disk galaxies with 
multi-hour, long-slit spectroscopy. The long-slit technique, 
although suffering from the much lower throughput in com- 
parison to the Fabry-Perot technique, allows us, for edge- 
on geometry, to take advantage of fields of view that cover 
the entire galaxy and greatly facilitates sky and continuum 
subtraction. We typically detect Ha emission to galactocen- 
tric radii of 1.5 R 25, and in certain cases up to 2 R25. In 
a previous paper (|Christlein fc Zaritskvl l2008h . we discuss 
the kinematic properties of these outer Ha disks and find 
the kinematics to be disk-like, with generally no indication 
of kinematic anomalies or higher velocity dispersions as one 
approaches the outer edge of the disk. A subsample of galax- 
ies with known optical warps has b een studied separately 
(|Christlein fc Bland-HawthorrJ [ioOSi ) to determine if kine- 
matic anomalies, such as breaks in the rotation curve, were 
associated with the onset of warps that could indicate on- 
going accretion processes of compact HI clouds or satellite 
galaxies as causes of the warp. We use data from both of 
these samples in this paper. 

We now turn our attention to the surface brightness 
profile of Ha emission in the outer disk. Our aim is to mea- 
sure the distribution of Ha emission in the outer disk and 
determine if the Ha emission profiles can be categorized as 
done for the stellar continua. We will determine whether 
there is indeed a break in the Ha surface brightness pro- 
file, and if so, whether this break can be associated with 
that in the stellar surface brightness. If there is a break, is 
it steep enough to constitute a truncation? And how is the 
stellar continuum, a measure of integrated star formation, 
distributed in comparison to the current star formation, as 
indicated by Ha? 

We describe our sample in ij2l In ^3.21 we discuss and 
compare the properties of Ha and stellar continuum emis- 
sion based on a composite of our 15 individual galaxy spec- 
tra, while in i]3.3l we discuss what constraints we can place 
on the profile shapes of individual galaxies. 
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2 DATA 

Our sample comes from three observing runs: a three-night 
run using FORSl at the VLT in November 2004, a three- 
night run with CMOS at Gemini-South in April 2005, and 
a three-night run using F0RS2 at the VLT in September 
2007. Target selection was, in all cases, for isolated, edge-on, 
late-type galaxies of several arcminutes in diameter, prefer- 
entially at a redshift of several thousand km/s, which places 
the Ha line in a relatively quiet window of the sky back- 
ground. An important difference among the runs is that the 
targets for the first two runs were chosen to be as mor- 
phologically undisturbed and regular as possible, whereas 
targets for the third run deliberately included several ob- 
jects with optical warps , most ly taken from the catalog by 
fSanchcz-Saavcdr a et al.l l|2003l ). At first, it may seem incon- 
gruous to include a set of warped galaxies in a long-slit ex- 
periment, because the disk material is eventually going to 
curve away from the slit at the onset of the warp. Ifow- 
ever, gi ven that nearly all galax:ie s might be warped to some 
extend jSanchez-Saavedra et al.l 12003 ). even if they appear 
"normal" in the available imaging, it is useful to include a 
subsample of known warped galaxies to provide the inter- 
nal control against which to view our results. Furthermore , 
warps typically set in around the R25 radius (Briggsl ll990l ). 
so that several important results from the undisturbed sam- 
ple within this radial distance may still be verified with the 
warped sample. 

All 15 targets are listed in Table [T] The Table lists the 
R25 radii adopted for this analysis, summarizes which tele- 
scope and instrument were used, and identifies whether the 
galaxy is k nown to be warped or not . Values of R25 are 
drawn from (|Lauberts fc Valentiin|[l989 h for all objects from 
the FORSl and GMOS-S runs a nd from the Third Reference 
Catalogue of Bright Galaxies (|de Vaucouleurs et al.l Il99ll ) 
for all objects from the FORS2 run. The projected optical 
diameter, 2R25, of all galaxies is at most half the slit length, 
so that each object is covered completely by the slit in a sin- 
gle observation, along with a substantial sky sample. In our 
analysis, we group the galaxies into an undisturbed and a 
warped sample. Two of the "unwarped" galaxies (ESO 340- 
G026, MCG -01-10-035) were observed during the Septem- 
ber 2007 run, which was primarily targeting warped objects. 

We describe the surface brightness profile as a function 
of projected radial distance from the galaxy center in units 
of the R25 radius so as to permit a comparison between 
galaxies of different physical sizes. Other characteristic size 
measures, such as the half-light radius, -R50, and the radius 
containing 90% of the total light, 7?9o, were also considered. 
However, all of these are gener ally strongly correlated (for a ll 
galaxies in our sample and in lLauberts &; ValentiinI (|l989l ). 
R25 ~R90 with a scatter of 15%, and with the exception of 
a single galax;y, R25 ~ 2XR50 to within 15%). We choose to 
use R25 because it is the only such quantity available in the 
literature for all of our galaxies and because it is a robust 
measure for the "edge" of the disk in this sample where 
inner-disk size measures can be heavily affected by dust. 



2.1 Data Reduction 

All raw spectra are processed using a standard IRAF 
pipeline. Flux calibration is carried out with spectropho- 



tometric standards taken each night; calibration factors for 
nights in a given run are consistent within 10%. Cosmic rays 
are identified by comparing a pixel value to the average of 
its neighbours within an annulus at distances between 3 and 
10 pixels and to the standard deviation in this annulus. The 
locus of pixels to be flagged as cosmic rays, depending on 
these quantities, is determined by visual inspection in a plot 
of the average neighbour pixel value versus the number of 
standard deviations. We subtract the sky background and 
continuum emission along the dispersion axis using the ap- 
propriate IRAF tasks. 

For each galaxy, we then extract kinematic data, using 
our own software. Line centroiding is carried out and veri- 
fied interactively for each individual row on the CCD. We 
fit and correct for the local Ha stellar absorption troughs, 
which are not removed by the sky subtraction or polyno- 
mial stellar continuum subtraction, with Gaussian functions 
(as now typically done f o r em ission line measurements, see 
iMoustakas fc KennicuttI (|2006i )). verify, and, if necessary, 
adjust the fit manually. We then transform all spectra onto a 
common coordinate frame, consisting of the projected sep- 
aration from the galaxy centroid in units of R25 and the 
wavelength offset from the interpolated wavelength corre- 
sponding to the mean Ha rotation velocity (i.e., we re- 
m ove the signature of ro t ation, for an example see Figure 8 
in^Christ lein Sz Zaritskvl (|2008 ^)). The transforming process 
does adversely affect the final signal-to-noise ratio, but has 
the advantage that all individual spectra can be compared 
on a consistent basis despite being taken with three different 
instruments and exhibiting different rotation curves. 

From the transformed spectrum, we measure the Ha 
flux within a tophat kernel that generously encompasses the 
entire line, including its deviations from the interpolated 
rotation curve. The flux is determined by integrating the 
counts across the entire kernel and renormalizing to units of 
erg s~^ cm~'^ arcsec"^. We determine uncertainties in the 
extracted flux by sliding the same tophat kernel along the 
wavelength axis in a "background" region without strong 
emission or absorption lines. We finally determine the stan- 
dard deviation of this residual flux over all background po- 
sitions and adopt that value as the uncertainty in the back- 
ground. We do not specifically include a term to account for 
Poisson uncertainties in the emission line flux, because our 
analysis focuses on a low-surface brightness regime where 
background errors are dominant. 

We calculate the stellar continuum brightness, using 
the same spectra without the stellar continuum subtraction, 
from the mean of background regions that lie adjacent to the 
Ha line and which are free of any strong emission or absorp- 
tion features. The uncertainty in this quantity is calculated 
as the error of the mean from the scatter of the individual 
pixel values. 



3 RESULTS 
3.1 Overview 

Before we begin our discussion, we set the stage by dis- 
cussing two key considerations. First, we need to define a 
set of models that we will test against the data. Our aim is 
to characterize the surface brightness profile of the Ha and 
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Table 1. Target list: R25 is the major axis radius of the 25 mag arcsec"^ isophote ((a): drawn from ESO-LV catalogue; (b): estimate 
based on acquisition image; close to RC3 value; (c): drawn from RC3). R"^ and R^^ are the outermost radii at which we detect Ha, 
according to the two definitions described in the text. 

stellar continuum emission quantitatively. For this purpose, The degree to which the composite represents a real galaxy 
we consider three basic hypotheses: rather than an amalgam of disparate objects will then be 

discussed when we examine individual galaxy profiles. 

• Hi: A broken exponential profile characterized by the 
equation 

/(r) = coe~^i''(l-e(r-ro))-hcoe~^i''«e"^2(r-ro)Q(^_^^) 3 2 Composite Hq and Continuum Surface 

Briglitness Profiles 

In this section, we focus on the general properties of the 
galaxies in our sample by superposing all individual spec- 
tra to create a composite spectrum. We take into account 
the fact that our sample consists of two subsamples, one 
of undisturbed, one of warped galaxies. In the case of the 
warped galaxies, we expect the galaxy isophotes to curve 
away from the slit position at and beyond the onset of the 
warp, which is likely to occur around a distance of R25 
(|Briggslll990t ) and should lead to a steep artificial drop of the 
surface brightness along the slit. Therefore, we will discuss 
the profile shapes of the two subsamples separately as well. 

In our default analysis, each galaxy is weighted equally, 
so that the resulting composite profile shape is effectively 
weighted by surface brightness. Although the average sur- 
face brightness scatters within only a factor of ~ 3 between 
most of the galaxies, there are outliers with very faint or 
bright surface brightness values, whose contributions to the 
composite will be affected by the choice of weighting; in 
particular, the measured surface brightness integral of IC 
2058 is almost three times as high as the second-brightest 
object. To test whether such outliers dominate our compos- 
ite profiles and thus bias our conclusions, we carry out an 
alternative analysis in which the contribution of each in- 
dividual galaxy is inversely weighted by the integral of its 
surface brightness profile over radius; in other words, the 
contributing spectra are normalized to a common mean sur- 
face brightness. 

In Fig. [T] we present the surface brightness profiles of the 
Ha and stellar continuum emission in our composite spec- 
tra. Even cursory visual inspection of the Ha and stellar 
continuum profiles shows that neither a single exponential 
nor a truncated exponential represent them well. Further- 
more, there are subtle differences between the undisturbed 



This is the most flexible model. It consists of two indepen- 
dent, but contiguous exponential profiles with a break radius 
of ro, characterized by a slope Ai in the inner, and A2 in the 
outer disk. The normalization is adjusted by the constant 
Co. The function 6 is the step function, i.e., 0{x) = 1 for 
X ^ 0, and otherwise. 

• H2'- A truncated exponential profile characterized by 

/(r) = coe-^i'-(l-e(r-ro)). (2) 

This surface brightness profile has a sharp cutoff at radius 
ro and is zero at larger radii. Free parameters are co, Ai, and 
ro. A sharp truncation is an idealized scenario, and even if 
it applied to the radial surface brightness profile, projection 
effects in a nearly edge-on sample would soften that edge. 
However, for typical parameters for slit width, disk scale 
height, scale length, and truncation radius, the effect on the 
profile shape is insignificant compared to the typical surface 
brightness fluctuations. 

• -ffa : A single exponential profile: 

fir) = coe-^i'-. (3) 

This model describes an unbroken exponential law, in accord 
with the Freeman Type I profile (|Freeman.l97Q ). Its only 
free parameters are cq and Ai. 

The second consideration is how we treat the data. We 
can either examine individual galaxy profiles in an effort to 
quantify the range of variation among profiles and whether 
any trends become evident, or combine the data into a single 
spectrum to enable a higher S/N analysis and, arguably, to 
reach general conclusions about the galaxy population. We 
choose to present the composite analysis first, to define what 
an "average" gala:xy might look like and address questions 
about the profiles from the highest possible S/N spectra. 
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R [R25] 

Figure 1. Ha and stellar continuum composite profile. The pan- 
els show, from top to bottom, the undisturbed, warped, and full 
sample. In each panel, the bold line represents the Ho surface 
brightness, the thin line the stellar continuum (surface brightness 
per A). The dotted line in the upper two panels shows the Ha 
profile from the bottom panel (full sample) to guide the eye. The 
dashed lines show the best fits to the Ha profiles. The shaded re- 
gion indicates the background error in the Ha flux. The vertical 
(flux) scale is linear up to y = 1, which corresponds to = 10~^* 
erg s~^ cm~^ arcsec"^; beyond, each integer step along the ab- 
scissa corresponds to a factor of 10. 

and warped composite profiles: The latter are slightly flat- 
ter overall and exhibit a particularly strong drop beyond 
~ I.I-R25, which may be a signature of the onset of the 
warp. For these reasons, we decide to fit all composite pro- 
files with broken exponential laws according to Eq. [T] It is 
possible that even a broken exponential law provides only 
an insufficient representation of the true profile shape, and 
that additional parameters would be required to model it ac- 
curately (e.g., a three-component broken exponential model 
for the warped objects); however, given the small size of the 
sample, we decide not to investigate more complex models 
in this paper. 

Given that the constituent galaxies may represent a 
range of profile shapes, the uncertainties in these fitting pa- 
rameters for the composite profile are likely to be dominated 
by scatter between the galaxies, rather than the measure- 
ment uncertainties. We therefore determine the uncertain- 
ties on these parameters by bootstrapping, i..e, by randomly 
resampling the set of profiles (two profiles per galaxy, rep- 
resenting the two sides) that we superpose to construct the 
composite, and fitting each of the resulting realizations of 
the composite spectrum separately. Because of the small size 
of the sample, we cannot guarantee that this procedure ade- 
quately samples the true variance in the parent population, 
but it provides a representation of the variance within the 
sample itself. 

We begin our analysis and comparison of the composite 
profile shapes by examining the break radius, ro, for both 



the Ha and stellar components. Fig.[2]shows the distribution 
of break radii recovered from our bootstrapping procedure 
for the undisturbed (solid contours), warped (dashed con- 
tours), and full (greyscale) samples. In the full sample, rep- 
resented by greyscales, the probability distribution for the 
Ha component exhibits a strong peak around ro ~ 0.7 and 
an extended tail towards larger radii, as well as a secondary 
peak for the stellar continuum break radius at ~ l.li?25. 
The 95% contours encompass the locus rna ~ r continuum', 
therefore, we cannot rule out that the break radii for the 
stellar continuum and Ha components are the same. The 
undisturbed and warped subsamples exhibit a qualitatively 
similar behaviour, with the difference that the secondary 
maximum for an Ha break radius around 1.1 R25 is much 
more pronounced in the warped than in the undisturbed or 
full samples. 

To further explore these results, we show the projec- 
tion of this probability distribution onto the r^™ axis in 
Fig. [3] All samples exhibit a strong peak in the probabil- 
ity distribution at ro ^ 0.7, as well as a secondary peak 
around ro ^ 1.0— 1.1. However, this secondary peak is much 
stronger in the warped than in the undisturbed sample; we 
therefore believe that it is related to the onset of the warp, 
where a drop in the surface brightness is expected. Never- 
theless, the evidence for the 0.7-R25 peak is clear even among 
the warped galaxies, thus justifying the decision to include 
both undisturbed and warped galaxies in our analysis. These 
results, modulo slight quantitative differences, are also ob- 
tained when weighting the individual spectra inversely by 
their integrated surface brightness in the composite (not 
shown in the figure, as the differences are small). 

We conclude 1) that the Ha break radius of 0.7 appears 
across both types of galaxies and so is not a result of warps, 
2) that this Ha break occurs in unwarped galaxies at r^°' — 
0-65lo;o6 (0-68lo.o2 for the full sample), and 3) that the 
break radii in the stellar continuum and Ha are consistent. 

The other important characteristic of the profiles that 
we consider in our comparison of Ha and stellar continuum 
surface brightness profiles are the inner and outer slopes. 
Fig.|4]shows the distribution of inner and outer slopes recov- 
ered from the bootstrapping procedure for both components 
(filled and empty circles for the Ha, crosses for the stellar 
continuum). Our discussion will focus on the full sample, 
but notable differences between the undisturbed and warped 
samples will be discussed. Our first conclusion from this Fig- 
ure is that the locus Ai = A2, indicated by a solid line in the 
bottom right corner of each panel, is inconsistent with both 
distributions; instead, A2 > Ai. This implies that the com- 
posite surface brightness profiles decline more steeply in the 
outer disk than in the inner one, i.e., both the stellar contin- 
uum and Ha profiles in the composite spectrum are of the 
Freeman Type II (the sub-exponential case) , and hypothesis 
Hi (Freeman Type I, the exponential profile) is ruled out. 

While the stellar continuum fits populate a fairly com- 
pact region in the Ai — A2 parameter space, fits to the Ha 
composite profile exhibit an extended tail. The maximum of 
the probability distribution for the full sample lies around 
Ai « 1.5 and A2 ~ 7, but a significant fraction of realiza- 
tions are fitted with much steeper slopes both in the inner 
and outer disk. Given our previous discussion of a dichotomy 
resulting from the combination of two samples with different 
selection criteria, it is natural to ask whether the two sub- 
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Figure 2. Probability distribution of the truncation radii tq 
of the Ha and stellar continuum components in the composite 
profile. Shaded regions show the distribution in the full sample, 
solid contours in the morphologyicall regular sample, and dotted 
contours in the warped sample. Contour lines are drawn at the 
p = 95% and p = 68% contours. Shaded areas represent parame- 
ter combinations with p > 5%, with darker shading representing 
the more preferred values. 



samples can be identified with distinct loci in Fig.|4l Indeed, 
selecting only those realizations of the bootstrapping proce- 
dure that yield best fits at ro < 0.9, i.e., those associated 
with what we consider the true truncation radius, produces 
fits solely in a well-confined region with Ai < 2.5 and A2 < 9, 
which we have indicated with filled circles. The remaining 
fits with ro > 0.9 are represented by empty circles. We will 
focus our analysis on the rg < 0.9 peak, which we believe 
to be uninfiuenced by the warps in several of our sample 
galaxies. 

A further important observation from Fig.|4]is that the 
outer disk slopes A2 for the Ha component are slightly dis- 
placed from those for the stellar continuum towards steeper 
slopes, indicating that Ha drops off slightly faster. Neither 
component is extremely large (> 10); with a mock catalogue, 
we have verified that a sharp truncation would yield values 
of A2 > 20. Projection effects, as noted earlier, will soften 
the outer edge of the surface brightness distribution only 
slightly and do not change this conclusion. Therefore, this 
result confirms that the truncation in neither component is 
particularly sudden, i.e., H2 is ruled out for the composite 
profile, as we had anticipated earlier based on visual inspec- 
tion alone. 

Finally, we turn our attention to the comparison be- 
tween the Ha and stellar continuum profiles. Irrespective of 
the infiuence of warped galaxies on the fit, we note that the 
two probability distributions are distinct from each other. 
In particular, fits to the stellar continuum are somewhat 
closer to the condition Ai = A2, meaning that the trunca- 
tion is relatively softer. The fits to the Ha component, on 
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Figure 3. Probability distribution for the break radius tq of the 
Ha surface brightness profile in the full sample (upper panel) and 
split up into the older sample of galaxies with undisturbed mor- 
phology (bottom panel, shaded histogram) and the newer sample 
of warped galaxies (bottom panel, unshaded histogram). The bi- 
modality of the probability distribution is much stronger in the 
warped sample; we therefore interpret the peak at ro = 1.1 as the 
onset of the warp. However, both samples independently show 
evidence for a break in the Ha surface brightness profile around 
ro = 0.7. 

the other hand, are marked by a significantly larger differ- 
ence between inner and outer slopes, i.e., the truncation is 
sharper. If we consider only the peak associated with the 
ro « 0.7 fits, we find the slope Ai of the Ha profile in the 
inner disk to be fairly shallow and much shallower than for 
the stellar continuum. In the outer disk, however, the situ- 
ation appears reversed: The values of A2 appear to indicate 
a slightly steeper slope for the Hq component than for the 
stellar continuum. 

The preceding observations also apply when we consider 
the samples of warped and unwarped galaxies separately; 
most differences are quantitative, but not qualitative. The 
most striking difference is the fact that the surface bright- 
ness profile of the warped galaxies is fitted with shallower 
inner slopes. One possible explanation is that the warped 
galaxies were selected to be highly edge-on systems, so that 
the surface brightness measured in the inner disk is strongly 
affected by dust in the galactic plane. However, since our 
analysis is not designed to probe the inner disks, we cannot 
pursue this observation further. More relevant to our inves- 
tigation, in the undisturbed subsample, the tail of the prob- 
ability distribution of fits to the stellar continuum profile 
extends to much steeper inner and outer slopes. This leads 
to a significant overlap between the probability distributions 
for the Hof and stellar continuum component. However, fur- 
ther investigation reveals that this tail is associated with 
fits with a very large break radius (> R25)', the apparent 
overlap in the fit parameters for the Ha (which all have a 
much smaller break radius) and stellar continuum profiles 
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Figure 4. Probability distribution for exponential slopes Ai and 
A2 in the inner and outer disk, respectively, for the Ha component 
(circles) and the stellar continuum (crosses). Panels show, from 
top to bottom, the undisturbed, warped, and full sample. Filled 
circles mark fits with Ha truncation radius ro < 0.9, and empty 
circles fits with ro > 0.9. The distribution for the Ha profile 
is bimodal, but the upper peak is associated with the truncation 
radius ro ~ 1.0 seen in the warped galaxy sample. The truncation 
is sharper in the Ho profile than in the stellar continuum. The 
line shows the locus Ai = A2, i.e., a single, unbroken exponential 
profile. 



is therefore caused only by the chosen projection, and the 
conclusion above that the two probability distributions are 
distinct from each other stiU holds. 

To better decide whether the steUar continuum and 
Ha components foUow significantly different spatial distri- 
butions or not, we examine an additional projection of the 
distribution of best-fit parameters. In Fig. [5] we examine 
the projection of the profile fit parameters onto the ro-A2 
plane, i.e., we plot the break radius versus the outer-disk 
profile slope. Again, filled dots represent realizations of the 
Ha composite profile, and crosses realizations of the stel- 
lar continuum profile. As previously, in order from top to 
bottom, the panels show the undisturbed, warped, and full 
samples. The fits to the Ha profile at ro « l.li?25 in the 
warped sample occur at very steep values of A2 > 15 and 
are thus outside the scale of this plot. In this figure, we 
clearly see that the overlap between the best-fit parameters 
for Ha and stellar continuum is negligible for all three sam- 
ples: For a given break radius ro, the Ha profile does not 
trace the stellar continuum, but declines significantly more 
steeply (but still with a finite slope). 

3.3 Individual Galaxies 

3.3.1 Statistical tests 

Throughout this section, we will test the likelihood of a given 
hypothesis describing the shape of the Ha or stellar contin- 



Figure 5. Probability distribution for the outer-disk exponential 
slope A2 vs. break radius ro for the Ha (dots) and stellar contin- 
uum surface brightness (crosses) profiles. Panels show, from top 
to bottom, the undisturbed, warped, and full sample. The two 
profile shapes are inconsistent with each other; for a given break 
radius ro, the Ha component drops more steeply in the outer 
disk. 



uum surface brightness profiles. These tests are based on 
quantifiers. However, it is not possible to quantify the abso- 
lute goodness of a fit using the estimator alone, because 
the scatter of the surface brightness measurements around 
the fit is dominated not by statistical measurement errors, 
but by intrinsic surface brightness fiuctuations of the emis- 
sion along the slit. These surface brightness fiuctuations are 
not refiected in our error bars; therefore, the formal quality 
of a x^ fit is always poor. Instead, we use a ratio 



(4) 



as our statistical quantifier, where X^{Hn) is the x^ value for 
a given hypothesis for the profile shape. The hypothesis Hi 
is our most fiexible one - it describes a broken exponential 
profile with a normalization factor co, a break radius ro, and 
slopes Ai for r < ro and A2 for r ^ ro. All other hypotheses 
Hn that we test are subsets of these four-parameter profiles, 
obtained by constraining one or more of the parameters to 
specific value. In general, the four-parameter model yields 
the best fit due to its having the largest number of degrees 
of freedom. This applies even in cases where Hn does pro- 
vide a correct description of the underlying surface bright- 
ness profile because small deviations from this shape due 
to noise and fluctuations are generally fitted better by the 
four-parameter hypothesis H\. However, a very large ratio 



indicates that the fit under hypothesis Hn is signifi- 



cantly worse than achievable with a four-parameter model, 
and that Hn has to be rejected. The virtue of this approach 
is that it is insensitive to the assumed uncertainties in the 
surface brightness measurement; the normalizations of the 
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uncertainties cancel out as we form a ratio of values. 
However, a complication lies in finding the critical thresh- 
old value for the ratio beyond which we can rule out that 
the improvement in the quality of the fit results simply from 
the ability to fit random fluctuations better with more free 
parameters, and is indeed indicative of a true preference of 
Hi over _ff„. 

To calibrate r^2 as a statistical measure of significance, 
we use a Monte Carlo simulation that generates mock sur- 
face brightness profiles conforming to the tested hypothesis 
Hn- The mock profile is created as a superposition of indi- 
vidual surface brightness fluctuations. The spatial shape of 
each profile is described by a Gaussian function with a fixed 
width. The number of fluctuations that are added up in each 
row of the simulated profile are determined from an assumed 
average amplitude and Poisson statistics. The actual flux of 
each Gaussian fluctuation is drawn from a Gaussian prob- 
ability distribution over logarithmic flux with a pre-defined 
width. In addition, we add fiuctuations with a fixed, small 
amplitude to simulate fluctuations in the background. The 
average amplitude, spatial width, width of the flux proba- 
bility distribution and background fluctuation amplitude are 
the parameters of this model; in addition, we have two pa- 
rameters controlling the size of the simulated uncertainties. 
We have run a series of simulations to find the parameters 
that recreate profile shapes most similar in appearance to 
the observed ones by fitting broken-exponential models to 
each mock profile and then using three statistics to compare 
them to the observations: 1) the ratio of the integrated resid- 
uals between 0.2 and 1.5 R25 to the total integrated profile 
in that range. 2) the cumulative residual between O.2R25 
and an outer radius x, normalized by the total integrated 
residual between 0.2 R25 and I.5R25. 3) the cumulative x^ 
between O.2R25 and an outer radius x, normalized by the 
total integrated between 0.2 R25 and I.5R25. The first 
of these statistics characterizes the magnitude of the fluc- 
tuations relative to the total flux. The second characterizes 
the radial behavior of the size of fluctuations. The third 
characterizes the relative goodness-of-fit as a function of ra- 
dius. Of particular importance is that, both in the mocks 
and in the real profiles, contributions to x^ come primarily 
from intermediate radii around r ~ 0.7, i.e., particu- 
larly sensitive to the region where we locate the break in the 
composite profile. 

We determine the ratio, , for a given hypothesis 
H„, and compare it to the scatter of r^2 in our mock pro- 
files. If r^2 is outside the range reproduced by 95% of the 
Monte Carlo realizations, we reject the hypothesis with 95% 
confidence. 

3.3.2 Individual galaxy Ha surface brightness profiles 

For our analysis of the Ha surface brightness profiles in the 
individual galaxies, we estimate the radial range of Ha de- 
tections in each galaxy. Although our fitting procedure in- 
cludes the entire radial range between 0.2 < R/R25 < 1.5, 
the position of the outermost detection additionally provides 
a lower limit on the position of the truncation in H2. 

To determine the radius of the outermost significant Ha 
detection, we integrate the total Ha flux between a given 
radius R and the largest radius considered in our analysis, 
1.5-R25. To determine the limiting radius Rum, we then use 




Figure 6. Ho (left) and stellar continuum (right) surface bright- 
ness profiles for individual galaxies, here from the morphologically 
undisturbed sample. Shaded regions in the Ha profile indicate 
background uncertainties. 

two alternative definitions. The first, Rfirm the radius of 
the outermost Ha detection with > 3cr significance between 
R and 1.5i?25. The second, Ru^n, is the smallest radius R 
for which the Ha flux is < 2a significant. The first defini- 
tion may include isolated Ha regions at large radial distance 
from the galaxy center, while the second definition measures 
only the extent of contiguous Ha emission and does not in- 
clude isolated sources in the outer disk. However, in practice, 
both definitions yield very similar results, except in the case 
of ESO 478- coil, where an Ha detection of low significance 
pushes the outermost radius to Rum ~ 1.45, while contigu- 
ous Ha emission is only detected out to Rum ~ 1.34. We 
list Riim for all galaxies in our sample in Table [1] 

We have not found obvious correlations between the ex- 
tent of the Ha emission and other global galaxy properties, 
such as the total Ha emission, apart from a very weak cor- 
relation between Rum and the integrated Ha flux beyond 
ro, which indicates that only galaxies with large Rum may 
have significant flux in this second, outer-disk exponential 
component. 

The Ha surface brightness profiles of the individual 
galaxies are shown in Fig. |B] (which consists of two parts) 
for the morphologically regular sample, and in Fig.[7]for the 
warped sample. By visual inspection, we identify a break in 
several of these galaxies, while others appear consistent with 
single exponentials. However, in general, there are large sur- 
face brightness fluctuations that preclude us reaching spe- 
cific conclusions from a cursory viewing. 

For this reason, a careful statistical analysis is required. 
We examine whether the individual galaxies are consistent 
with specific hypotheses for the shape of the composite pro- 
file in the manner described in i i3.3.1l We calculate the 
ratio between fits with the specific hypothesis and a more 
general four-parameter, broken exponential fit, and then de- 
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Figure 6. (continued) 
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Figure 7. Ho (left) and stellar continuum (right) surface bright- 
ness profiles for individual galaxies, here from the warped sample. 
Shaded regions in the Ha profile indicate background uncertain- 
ties. 



termine, using a Monte Carlo simulation, if this ratio indi- 
cates that a four-parameter model provides a significantly 
better fit and that the tested hypothesis must be discarded. 
In all cases, we construct the input profile shape for the 
Monte Carlo simulation on the basis of the best-fit param- 
eters for the composite profile under the given hypothesis 
(broken, truncated, single-exponential). We find no qualita- 
tive global differences when using the two different weighting 
schemes in our construction of the composite spectra, al- 



Name 


H2 


Hi 


Hla 


H4 


H2a 


Warp? 


ESO 201- G 022 


0.00 


0.00 


0.06 


0.01 


0.00 


N 


F.SO 2QQ- n 018 


0.04 


0.69 


0.91 


0.19 


0.02 


N 


F.SO '^2S- a O'^S 


0.00 


0.00 


0.00 


0.00 


0.00 


N 


HjOkJ OOU Vjj \j^o 


0.02 


0.00 


0.44 


0.00 


0.00 


N 


ESO 385- G 008 


0.03 


0.08 


0.82 


0.00 


0.00 


N 


ESO 478- G Oil 


0.06 


0.92 


0.90 


0.05 


0.04 


N 


ESO 340- G 026 


0.00 


0.03 


0.80 


0.65 


0.00 


N 


IC 2058 


0.07 


0.37 


0.67 


0.00 


0.00 


N 


IC 4393 


0.02 


0.01 


0.06 


0.00 


0.00 


N 


MCG -01-10-035 


0.02 


0.68 


0.77 


0.12 


0.02 


N 


ESO 184- G 063 


0.02 


0.00 


0.39 


0.00 


0.00 


Y 


ESO 473- G 025 


0.03 


0.00 


0.24 


0.01 


0.01 


Y 


NGC 259 


0.03 


0.86 


0.89 


0.01 


0.03 


Y 


UGC 12423 


0.01 


0.21 


0.47 


0.04 


0.00 


Y 


UGCA 23 


0.00 


0.00 


0.89 


0.11 


0.00 


Y 



Table 2. Fraction of simulations that exceed the observed 
Xjjn/x\i ratio for a given galaxy, where Hn is one of the fol- 
lowing: H2'- trunated exponential, Hy. single exponential, Hia'- 
broken exponential with tq = 0.7, H^ : scaled version of composite 
best fit, H2a'- truncated exponential extending to last significant 
Ha detection. 



though a few individual galaxies shift from being consistent 
with a given hypothesis to being inconsistent and vice versa. 
For the remainder of this discussion, we use the composite 
constructed from the unweighted combination of individual 
spectra. 

Using the lessons learned from the composite spectrum 
and the details evident in the individual spectra, we mod- 
ify our set of models to include a version of Hi, the broken 
exponential, that sets the break radius to 0.7 (model Hia), 
a version of H2, the truncated exponential, that sets the 
radius of truncation to be that of the outermost Ha detec- 
tion (model H2a), a single-exponential model {H3), and a 
new model that uses the best-fit composite profile (model 
H4 with free parameter co to allow renormalization) . The 
fit quality is evaluated only between 0.2 and 1.5 -R25. 

We present the results of this analysis in Table (2] The 
entries in the Table give the fraction of realizations for which 
the ratio Xhu/Xhi exceeds the value obtained from the ob- 
served profile. Therefore, low values indicate that hypothesis 
Hi provides a significantly better fit to the surface bright- 
ness profile than the alternative. We highlight entries < 0.05, 
i.e., those that formally indicate that the considered hypoth- 
esis is inconsistent with the observed profile shape at the 2a 
level. 

We now discuss the results in terms of specific questions. 

3.3.3 Are the individual galaxies consistent with a 

truncated exponential surface brightness profiles? 

Hypothesis H2 is that of an exponential profile with a sharp 
truncation. Most galaxies in our sample, whether warped or 
not, are inconsistent with the truncated profile shape. The 
two cases that are marginally consistent with a sharp trun- 
cation are ESO 478-COll and IC 2058. Visual inspection 
supports this for IC 2058; in the case of ESO 478-GOll, 
the best-fit truncation radius is at the very large radius of 
1.347?25, so that the truncated profile is essentially indistin- 
guishable from a single exponential shape. 
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We now restrict this hypothesis. In many of our sample 
galaxies, Ha is observed to large radii, as Table [T] shows. 
Such detections are formally inconsistent with an exponen- 
tial profile truncated at smaller radii. Therefore, we impose 
the additional constraint that the truncation radius be lo- 
cated just beyond the radius of the outermost Ha detection 
and choose Rum for this purpose. In Table (2] the results of 
this analysis are shown in column H2a - None of the galaxies 
in our sample are described well by a profile truncated at 
the last significant Ha detection. 

Lastly, perhaps a sharply truncated profile is ruled out 
by this analysis only because the truncation is softened in 
projection due to the nearly-edge-on configuration. Simple 
simulations of the projected profiles that would be obtained 
from such a configuration show that, even for an almost 
perfectly edge-on geometry, the softening effect towards the 
edge of the Ha disk would be fairly subtle and unlikely to 
exceed the surface brightness fluctuations observed in the 
real spectra. Therefore, this analysis conflrms that a broken- 
exponential profile with a finite outer-disk slope provides a 
significantly better fit to the individual galaxy profiles than 
a sharply truncated single exponential law. 

3.3.4 the individual galaxies consistent with an 

unbroken exponential profile? 

Next, we consider H^, the hypothesis that the Ha sur- 
face brightness profile is described by a single, unbro- 
ken exponential. This model corresponds to the Freeman 
Type I profile, which is obs erved to hold for the stellar 
continuum of some gal a xies dBland-Hawthorn et al.l l2005l : 
IPohlen fc Truiillol I2OO6I : lErwin. Pohlen fc BeckmanI 120081 ) 
out to extremely large radii. 

Our sample is dominated by galaxies with broken ex- 
ponential profiles in the stellar continuum (Figures |6] to [T} , 
where visual inspection identifies only few objects as being 
consistent with a single exponential profile. For about half 
of the Ha profiles in our sample, whether they represent 
normal or warped galaxies, we also rule out the unbroken 
exponential profile shape (i/a), as Table [2] shows. The ex- 
ceptions among the unwarped galaxies may be of particular 
interest because the existence of a critical threshold surface 
gas density beyond which the specific star formation rate in 
a galaxy changes abruptly or even dr ops to zero has been 
sugge sted repeatedly in the literature ijMartin fc KennicuttI 
l200lh . The objects for which we cannot conclusively rule 
out an unbroken exponential law are ESO 299-G018, ESO 
385-G008, ESO 478-GOll, IC2058, MCG -01-10-035, NGC 
259, and UGC 12423. In some cases, this can be attributed 
to the strong fluctuations, which lower the quality of even 
the broken-exponential flt and reduce our ability to discrim- 
inate between different hypotheses. Visual inspection does 
single out two of these objects, ESO 478-GOll and NGC 259, 
as the best candidates for an unbroken exponential profile; 
these are also the two objects best fitted under H3,. This is 
surprising for NGC 259, which is listed as a warped galaxy; 
however, within our warped sample, this object is the least- 
inclined, and the classification as a warped object may be 
due to a prominent spiral arm, rather than to an actual verti- 
cal bending of the disk, which may explain why we can trace 
a regular, possibly unbroken Ha profile out to ^ 1.35 x R25. 

We conclude that, just like the composite profile, the in- 



dividual galaxy Ha profiles are generally inconsistent with 
unbroken exponentials out to large radii, but that a signifi- 
cant number of exceptions might exist in our sample. 

3.3.5 Are the individual galaxies consistent with a 
universal break radius? 

The third question we address is whether the individual 
galaxies might all share a common break radius, ro « 0.7, 
as measured in the composite profile. Is this radius a uni- 
versal property of late-type galaxies, or simply an average 
for a population with an intrinsically large scatter? We ad- 
dress this question by starting with Hi, and then fixing the 
break radius at ro = 0.7, but leaving cq, Ai, and A2 as free 
parameters. We refer to this model as Hia. The result, pre- 
sented in the Hia column of Tabled shows that all but one 
(ESO 323-G033) of our galaxy profiles are consistent (< 2a 
discrepant) with broken exponential profiles with a break 
radius of ro = 0.7. ESO 323-G033 bears the distinction of 
being consistent with neither of the four alternative profile 
shapes we are considering here. Because a single inconsis- 
tency at the 2a level in a sample of 15 is within statistical 
expectations, we conclude that the hypothesis of a "univer- 
sal" break radius for late-type galaxies like we have selected 
is consistent with these data. 

3.3.6 Are the individual galaxies consistent with the 
composite profile? 

Extending the previous analysis, we now examine whether 
the composite profile calculated in H3.2I is itself an accept- 
able representation of all individual galaxies in the sample. 
Again, we create mock profiles, this time with the composite 
spectrum for which we allow a renormalization. The results 
(Table [2I1 show that the Ha profiles of most of our galaxies 
are inconsistent with the composite; there are only four ex- 
ceptions: one from the undisturbed sample (ESO 299-G018) 
and three from the warped sample (ESO 340-G026, MCG 
-01-10-035, UGCA 23). We conclude that there are real vari- 
ations among the galaxies and that the composite spectrum 
is not necessarily a good representation of a galaxy cho- 
sen at random, with the exception of the break radius at 
ro = 0.7i?25. 



4 DISCUSSION 

We find that more than half of our galaxies do not have 
a Freeman Type I Ha surface brightness profile (a single, 
unbroken exponential). Visual inspection demonstrates that 
the profiles are rather of the downbending (sub-exponential) 
Type II class, and are thus qualitatively following the be- 
haviour of the stellar continuum, i.e., their surface brightness 
profiles may be described as broken exponential profiles with 
a characteristic break radius and a steeper slope in the outer 
disk than in the inner disk. The Ha break radius at ro ~ 0.7 
R25 appears surprisingly well-constrained and is observable 
both in the unperturbed and the warped subsamples. 

Observations of a break in the Ha surface brightness 
profile are not new. Breaks or trunca tions in the Ha pro- 
file have been noted for many years (cf. iMartin fc KennicuttI 
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I2001I ). and a star formation threshold, i.e., a critical gas sur- 
face sensitivity below which star formation cannot proceed 
or proceed only very inefficiently has often been suggested as 
the most likely explanation. Other breaks or truncations are 
expected at very large radii, where ionization by the cosmic 
ultraviolet background may become a significant contributor 
to outer-disk Ha emission and where the baryonic disk ends; 
however, both of these effects are likely occur at or beyond 
1.5 -R25, which is a typical outer radius of HI disks, and at 
surface brightness levels fainter than what we can reliably 
probe with these data. What is new here is that we quantify 
the radial behavior of the emission beyond this break radius 
and out to large radii, ~ 1.5-R25. 

It is tempting to identify the clear break in Ha with 
the Ha, and possibl y star formation, threshold well-known 
from the literature (|Martin fc Kennicut3l200ll '). but if so, 
then this "threshold" is clearly not the end of star forma- 
tion because Ha emission continues with some regularity 
beyond. Furthermore, if the threshold position depends on 
local physical conditions, it is remarkable how well-defined 
the location of ro is with respect to f?25. 

It is also tempting to identify the break radius in the 
stellar continuum t hat we observe with that observed by 
iPohlen et all (|2002l ). and then further associate it with our 
observed Ha break radius because the two are consistent 
within 2cr in our data. However, the outer-disk exponential 
slopes are different for the two components, with the stellar 
continuum being more extended, suggesting that there is not 
a 1:1 correspondence between the surface brightness profiles 
of the two components. The physics determining the two are 
therefore more complex than a simple relationship between 
star formation and stellar populations. 

If Ha is a linear tracer of star formation, the cur- 
rent stellar continuum surface brightness distribution can- 
not have been produced by in-situ star formation, unless star 
formation was more extended in the past. Because the cur- 
rently favored paradigm is that galaxies form "inside-out" 
(|Truiillo fc Pohlenll2005l : I Williams et "al]|2009l ). we consider 
it more likely that either the outer-disk stellar population 
has resulted from a redistribution of stellar mass (for ex- 
am ple, as a result of tidal interactions with other galax- 
ies dKazantzidis et al.l[2008l ) or secular dynamical evolution 
l|Roskar et al.ll2008l )). or that Ha is an increasingly biased 
tracer of star formation towards large radii. 

No less important than the observation of a break in 
the Ha surface brightness profile between the inner and the 
outer disk is the complementary observation that Ha emis- 
sion is not sharply truncated. We observe Ha emission as 
far out as 1.5ii25, and find its surface brightness profile to 
continue with some regularity, possibly even following an ex- 
ponential law, but with a steeper slope. The observation of 
such extended emission is consistent with observations in re- 
cent years of extended UV emission around seve ral galaxies 
iThilker et al.ll2005l : [Zaritskv fc Christleinll2007l ). It demon- 
strates that star formation can occur even in the extreme 
outer-disk regime, where gas surface densities are presum- 
ably very low. At the same time, the presence of the break 
in the Ha profile suggests that, even if star formation is 
not completely suppressed in the outer disk, it appears to 
proceed under different physical conditions. 

The finding that the outer disk UV emission, which is 
also an indicator of recent star formation, does not exhibit 



a break (|Boissieij|2007al lbh. has given rise to the suggestion 
that the star formation rate profile is also unbroken. A pos- 
sible explanation of the divergent behaviour between the Ha 
and UV profiles is that the truncation of Ha emission is a 
stochastic effect due to the number of star formation regions 
with sufficiently massive, ionizing stars dropping below unity 
(|Boissieill20"07al '). However, in our observations, which reach 
much higher sensitivies than narrow-band imaging studies 
and superpose multiple galaxies, we detect emission at ex- 
tremely high levels of confidence, at and beyond the hypo- 
thetical truncation radius. Therefore, there are numerous 
ionizing sources in the outer disk. Nonetheless, we see a 
clearly-defined break in the Ha surface brightness profile 
both in individual galaxies and in the composite; we con- 
clude that the break in the Ha profile is real, but that it 
does not mark the end of star formation or the end of the 
formation of sufficiently massive stars. 

The apparent disagreement with the UV surface bright- 
ness profile is thus puzzling. If we interpret the UV observa- 
tions to mean that star formation is not subject to a sudden 
break, then Ha must be a biased tracer of star formation 
and suppressed at large radii. The ratio of UV to Ha emis- 
sion must change continuously with radius. One potential 
cause of such an effect is an increasingly bottom-heavy ini- 
tial mass function with radius, either as a result of a real 
change in the functional form or from statistical effects re- 
lated to the increasingly poorer sampling of the upper stellar 
masses as total mass decreases. A scenario in which the in- 
tegrated galactic IMF (which is probed here) may vary with 
environment as a result of variations of the mass function 
of molecular clouds across the galaxy, while the shape of 
the IMF itself would be preserved, has been presented by 
iKroupal (120071 1. 

Another potential cause is a transition from ionization- 
bounded to density-bounded HII regions. If a larger fraction 
of ionizing photons escape without being converted to Ha at 
larger radii the UV to Ha ratio would change. Within inner 
galaxy disks the effect seems to go the other way, with high- 
luminosity star-forming regions, rather than low-lum inosity 
ones, appearing to be density bounded (Bcckma n et al.l 
[2000), but the situation in the outer disk could be quite 
different. 

A final explanation for the increasing bias between 
stellar continuum and Ha emission in the outer disk is 
that outer-disk star formation may be an intermittent phe- 
nomenon, and short, but correlated, episodes of star forma- 
tion in the past may have created the present extended stel- 
lar population before returning to the low-level state that we 
observe in our sample. A similar scenario, involving bursts of 
star formation across the entire disk, has been proposed by 
iBoissieJ (|2008t ) based on GALEX observations of UV colors 
of outer galaxy disks. However, this scenario does not ap- 
pear to readily explain why the UV surface brightness pro- 
file would not exhibit a truncation, unless the galajcy was 
observed shortly after such a burst, so that the UV outer 
disk would still be dominated by stars formed during the re- 
cent star formation episode, while Ha emission, which traces 
star formation only over very short timescales, has already 
aba ted. 

iBakos. Truiillo fc Pohleiil (|2008l 'l find that the stellar 
surface mass density also exhibits a much less pronounced 
break than the optical light of the stellar continuum. They 
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suggest that a change in the composition of the steUar popu- 
lation, rather than an actual break in the mass distribution, 
is responsible for the observed break in the optical contin- 
uum surface brightness. In this context, it is surprising that 
both the stellar surface mass density and the integrated, 
recent star formation rate (as measured by UV emission) 
supposedly show no or only weak breaks, but both stellar 
continuum and Ha do. 

As exceptions to the rule, those galaxies where we find 
no evidence for a truncation or break in the Ha emission 
deserve particular attention, as any theory attempting to 
explain the presence of a break in the Ha profile — for 
example, by invoking a critical gas surface density threshold 
below which star formation is inefficient — must also be able 
to provide an explanation for objects that apparently violate 
the rule. Of the objects that are formally consistent with an 
unbroken exponential profile, NGC 259 and ESO 478-GOll 
may be interesting cases for follow-up studies. 

A final consideration must be given to the effect of dust 
on the Hq surface brightness profiles, and, in particular, to 
the question whether dust absorption could influence our 
measurement of the characteristic break radius and possi- 
bly even our classification of the composite Hq profile as 
a (sub-exponential) Type H, rather than a (single exponen- 
tial) Type I. If dust extinction sets in around this radius and 
increases continuously towards the galaxy center, absorbing 
an increasing fraction of Ha flux, it may artificially flatten 
the profile slope and possibly even create the false impres- 
sion of a broken exponential shape. This scenario, however 
would imply that the unabsorbed Hq profile could be recre- 
ated by extrapolating the outer-disk Ha profile back towards 
the center, implying that the real Ha surface brightness in 
the inner regions would be more than an order of magnitude 
higher than observed here. As dust shielding is unlikely to be 
uniform (even in cases of nearly perfect edge-on alignment, 
the slit intercepts regions both on and off the central dust 
lane), at least some sight lines should permit glimpses of 
these substantially higher surface brightness levels, leading 
to a larger dynamic range in the fluctuations of the inner- 
disk surface brightness profile. This is not the case; even 
the highest Ha peaks remain far below the levels expected 
from interpolating the outer-disk profile, confirming that the 
break is indeed real. Dust shielding may, however, plausibly 
modify the measured value of the inner-disk profile slope by 
smaller amounts. 

An additional argument against dust obscuring a large 
fraction of the Ha fiux comes from the rotation curves. 
If dust played a significant role up to the break radius of 
0.7 R25, contributions to the observed Hq flux would be 
weighted towards gas on the outskirts of the galaxy along 
the edge closer to the observer. We would therefore expect 
the rotation curve to rise very slowly in the inner disk and 
not flatten before it reaches this point. However, most rota- 
tion curves in this sample are already fiat or close to fiat 
well before reaching 0.7 R 25 fehris tlein fc Zaritskvl l2008l : 
IChristlein fc Bland-Hawthorn 20081. indicating that the Ha 
emission at this point already samples the kinematics of the 
gas along the entire sight line, and that the observed kine- 
matics and fiux cannot be strongly biased by internal dust 
extinction. 



5 CONCLUSIONS 

We have presented a study of Ha emission and stellar con- 
tinuum in 15 low-redshift, late- type, edge-on galaxies. The 
use of deep long-slit spectroscopy, rather than narrow-band 
imaging, allows us to probe Ha emission to very faint levels, 
^-|^Q-i8 gj.g g-i arcsec"^. Ha emission is traced out 

to 50% beyond the R25 radius, roughly twice as far as typ- 
ical Ha rotation curves in the literature. It is thus possible 
to probe the outer galactic disks, which have hitherto pri- 
marily been a domain of radio astronomy, in Ha emission, 
and thus simultaneously gain information about kinemat- 
ics, star formation, metallicity, and stellar continuum with 
arcsecond-scale spatial resolution. 

In the present paper, we have focused on the properties 
of the surface brightness profile of Ha emission and com- 
pared them to those of the stellar continuum. Our conclu- 
sions are: 

• In this sample, which is dominated by Freeman Type 
II ("sub-exponential") profiles in the stellar continuum (a 
broken exponential profile with a steeper slope in the outer 
than the inner disk), Hq emission in the composite profile 
also follows a Type II profile. For the Ha composite profile, 
we can rule out both an unbroken single exponential profile 
as well as a truncated profile, for which the fiux drops to 
zero at a given break radius. 

• There is a well-defined break radius in the composite 
Ha surface brightness profile at ro ~ 0.7 R25. 

• The Hq and stellar continuum distributions may be 
consistent with a single break radius ro, but the stellar con- 
tinuum profile drops more slowly in the outer disk. 

• The presence of a clearly defined break in the Hq sur- 
face brightness profile, despite the fact that our observations 
of the outer disk slope are clearly not limited by Poisson 
noise, indicates that any apparent deficit of Hq emission in 
the outer disk is not, in general, a stochastic effect arising 
from low number statistics. The break in Hq is real, but it 
is not a truncation, as evidenced by the common presence 
of low- level outer-disk Ha emission to ^ 1.5ro. 

• For most objects in our sample, a truncated exponential 
as well as an unbroken exponential can be ruled out even on 
an individual basis, thereby requiring a broken exponential 
profile. All but one of our galaxies are individually consis- 
tent with a broken exponential profile, with a fixed Hq break 
radius at ro — 0.7ii25, showing that the conclusions drawn 
from the composite profile describe not just a meaningless 
average, but are representative of at least a significant frac- 
tion of individual galaxies. 
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